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Abstract 

The study of impurities in low dimensional antiferromagnets has been a very active field in magnetism ever since 
the discovery of high temperature superconductivity. One of the most dramatic effects is the appearance of large 
Knight shifts in a long range around non-magnetic impurities in an antiferromagnetic background. The dependence 
of the Knight shifts on distance and temperature visualizes the correlations in the system. In this work we consider 
the Knight shifts around a single vacancy in the one and two dimensional Heisenberg model. 
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Doping and impurity effects in low-dimensional 
quantum antiferromagnets remain of strong interest 
in the condensed matter physics community, spurred 
by high-Tc superconductivity and other exotic efl^ects 
in those systems. One of the most basic questions to 
ask is how an antiferromagnetic systems responds to a 
uniform applied magnetic field and how this response 
is changed in the presence of impurities. This question 
will be answered in detail for a two-dimensional (2D) 
spin- 1/2 Heisenberg model. The results will then be 
compared to the one-dimensional spin-chain model 
which does not order at low temperatures. 

Let us first consider a generic antiferromagnet, i.e. a 
collection of spins of size S which are assumed to or- 
der antiferromagnetically at low temperatures. This 
implies two sublattices A and B with long range cor- 
relations throughout the lattice at low temperatures. 
The magnetic moments are correlated in parallel on the 
same sublattice, but antiparallel to the other sublat- 
tice. The size of the order parameter may be reduced 
by temperature or quantum fluctuations, but it is as- 
sumed to be non-zero. The effect of an applied mag- 
netic field B can be intuitively understood as described 
by standard textbooks [1]. If the field is perpendicular 
to the order (transverse field), all spins can tilt slightly 
towards the field as shown in Fig. 1. A uniform mag- 
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Fig. 1. The effect of a magnetic field B on a generic antiferro- 
magnet. 

netization is induced in the entire sample with a finite 
magnetic susceptibility x = dm/dB — —d^F/dB^ at 
low temperatures. Note, that if the field was applied 
parallel to the antiferromagnetic order, the response 
would be very small, since it costs more energy to in- 
duce a longitudinal magnetization. In an isotropic an- 
tiferromagnet, the order is therefore always aligned 
transverse to the field as shown in Fig. 1. At non-zero 
temperatures, spin-waves become excited, which can 
be polarized with the field. Therefore the susceptibility 
increases with increasing temperature in the ordered 
phase. On the other hand, at very high temperatures 
in the non-ordered phase the susceptibility is well de- 
scribed by the Curie- Weiss-law S{S + l)/3(r + 6), 
i.e. decreasing with temperature. The typical suscep- 
tibility of an antiferromagnet therefore shows a broad 
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Fig. 2. Generic magnetic susceptibility for an antifcrromagnct. 

correlation maximum as shown in Fig. 2. Even in anti- 
ferromagnetic models which do not order at low tem- 
peratures this correlation maximum is well established, 
as for example in the spin-1/2 chain[2]. In this system 
the entangled quantum state gives rise to exotic effects 
such as a diverging derivative of the susceptibility with 
respect to Tas T^ 0[3]. 

We will now considers what happens to the indi- 
vidual spins when a vacancy is introduced into the 
antiferromagnetic order on sublattice B as shown in 
Fig. 3. Because of the antiferromagnetic order all spins 
on sublattice A form one large magnetic moment of 
size SN/2, while all spins on sublattice B amount to 
a magnetic moment of S{N/2 — 1). The total system 
has therefore an effective net classical moment of size 
S which can align with the magnetic field, following 
a Curie-law with diverging susceptibility S'^ /3T. This 
means that all spins on the sublattice A have a diverg- 
ing susceptibility of S^ /ST and all spins on sublattice 
B have a negative diverging susceptibility of —S'^/3T. 
This alternating Curie-type response is much larger 
than the uniform canting of the spins at low tempera- 
tures discussed above. One single impurity is therefore 
sufficient to induce an antiferromagnetic longitudinal 
order throughout the lattice in the presence of a mag- 
netic field. So far we have neglected that in a quantum 
antiferromagnet the order parameter m is reduced by 
quantum fluctuations (m = 1 represents the maximum 
alignment of all spins on the respective sublattices) . 
Since the induced longitudinal order stems from the 
existing transverse order in the system, it must also be 
proportional to mS^ /ST (m ~ 80% for the spin-1/2 
Heisenberg AF on a 3D cubic lattice, and m ~ 60% 
on the 2D square lattice). Therefore, we can write ap- 
proximately for the local response around an impurity, 
according to this intuitive picture 



X(r) = Xo + (-1) 



r^+ry + l 



ST' 



(1) 
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Fig. 3. The effect of a magnetic field B on an antiferromagnet 
with one single impurity (vacancy). 

where xo is the susceptibility per site of the pure 
system. The induced response is alternating (Fig. 3) 
while the uniform response remains largely unchanged 
(Fig. 1). 

We will now test this picture for the 2D spin-1/2 
Heisenberg model, 



H = J ^ S, 



(2) 



(i,j) 



where (i, j) denotes nearest neighbor sites on a periodic 
square lattice. This model is known to exhibit antiferro- 
magnetic order as T ^ 0. However, for any flnite tem- 
perature the order is destroyed by quantum fluctua- 
tions. Among the ordered antiferromagnets, this model 
is therefore at the extreme borderline to a quantum en- 
tangled state. The quantum Monte Carlo program we 
developed uses the loop algorithm in a single cluster 
variety implemented in continuous time[4,5,6,7], which 
gives efficient and fast updates even at very low tem- 
peratures. 

According to Eq. (1) the local response around a 
static vacancy 



x{r) = pj2(^:s:} 



(3) 



can be separated into a sum of uniform and staggered 
parts on the lattice 



x(r)=Xuni(r) + (-ir+'-''x.tag(r), 



(4) 



the amplitudes of which are both slowly varying on the 
scale of one lattice spacing. In order to extract those 
two components we numerically extrapolate the data 
on the even sublattice to the odd sublattice and vice 
versa and define 



Xuni(r) = 
Xstag(r) = 



Xcvcn(r) +Xodd(r) 

2 
Xcvcn(r) -Xodd(r) 



(5) 
(6) 



The results for the staggered and uniform parts are 
shown in Figs. 4 and 5 for T — 0.05J. The uniform 
part drops off very fast to the limiting value xo , but is 
strongly enhanced around the impurity. The staggered 





Fig. 6. The ratio 12TXoo as a function of /3 = l/T. 



Fig. 4. (Color online) Xstag(r) for T=0.05J 




Fig. 5. (Color online) Xuni (r) for T=0.05J 

part also approaches a limiting value Xoa ~ 0.6|y, 
which has the expected Curie-behavior as shown in 
Fig. 6 in agreement with m ~ 0.6 for the 2D Heisen- 
berg model[8,9]. However, a broad peak around the 
impurity remains, which appears to be largely tem- 
perature independent. Therefore Eq. (1) appears to 
be indeed valid for longer distances from the im- 
purity, while for shorter distances quantum effects 
dominate. We expect that this picture becomes more 
and more accurate for larger spin models in higher 
dimensions. Our simulations for the 2D model are in 
full agreement with previous calculations, where the 
impurity susceptibility of the entire system has been 
considered[10,ll,12,13,14,15,16,17]. For more details 
on the effect of more than one impurity, see Ref. [18]. 
We now wish to compare the situation to a model 
which does not order at low temperatures, such as 
the ID spin-1/2 Heisenberg model. In this case, va- 
cancies cut the chain, so that we need to consider a 
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Fig. 7. The local response for a semi infinite chain at T — 0.05 J. 

semi-infinite chain. This problem has been considered 
before[19], with some surprising results. A large alter- 
nating response is also induced by the impurity, but 
this increases with the distance from the edge. Finite 
temperatures, finite fields, or finite system sizes will 
limit the range of the alternating part, but generally 
the maximum alternating response is not closest to the 
impurity site. A typical response is shown in Fig. 7 for 
T = 0.05. It is clear that such a complicated pattern 
is an indication of a collective state in this quantum 
many body system, which gives some indication of the 
nature of the valence bond state [20]. 

In conclusion we have analyzed the local response to 
a uniform magnetic field in low-dimensional spin-1/2 
antiferromagnets with one vacancy. For the 2D model 
an intuitive picture of a long range antiferromagnetic 
order describes the results well far away from the impu- 
rity site. Only the more local enhancement is specific to 
the model and must be attributed to quantum effects 
or corrections from spin- wave theory. We therefore can 
argue that the local response is always accurately de- 
scribed by Eq. (1) for larger distances from a vacancy 
in any ordered antiferromagnet, where m is the order 



parameter, which typically depends on spin, dimension 
and temperature. Therefore a single impurity induces a 
large Curie-divergent alternating response throughout 
an ordered antiferromagnet, if a field is applied trans- 
verse to the underlying magnetic order (this is auto- 
matically the case in an isotropic model) . The response 
in the immediate vicinity is specific to the model, but 
can sometimes also be intuitively understood, e.g. by a 
valence bond basis[20]. The situation is quite different 
for the ID quantum antiferromagnet, where strongly 
entangled quantum states dominate the picture. 

Our results will have direct consequences on NMR 
and /iSR experiments on doped antiferromagnets[21,22,23 
For the one-dimensional case the exotic boundary ef- 
fects have already been confirmed, which are manifest 
through NMR satellites with a characteristic l/\/r 
dependence [24]. For ordered antiferromagnets it is im- 
portant to distinguish between the ordered sublattice 
magnetization [25] and the field induced staggered 
magnetization, which we have described here. The 
field induced effects should show a more dramatic 
Curie-like temperature dependence and are assumed 
to be perpendicular to the sublattice magnetization. 



[14] S. Sachdcv, J. Stat. Phys. 115, 47 (2004) 

[15] M. Vojta, C. Buragohain, S. Sachdcv, Phys. Rev. B 61, 
15152 (2000) 

[16] O.P. Sushkov, Phys. Rev. B 68, 094426 (2003) 

[17] S. Sachdcv, M. Vojta,, Phys. Rev. B 68, 064419 (2003) 

[18] F. Anfuso and S. Eggcrt, preprint iCond-mat/051 10011 
(2005). 

[19] S. Eggert and I. Affleck, Phys. Rev. Lett. 75, 934 (1995). 

[20] G.B. Martins, M. Laukamp, J. Riera, E. Dagotto, 
Phys. Rev. Lett. 78, 3563 (1997); M. Laukamp, 
G.B. Martins, et al, Phys. Rev. B 57, 10755 (1998) 

24l[21] A.V. Mahajan,et al., Phys. Rev. Lett. 72, 3100 (1994). 

[22] S.T. Ting, et al. Phys. Rev. B 46, 11772 (1992). 

[23] M. Corti, et al. Phys. Rev. B 52, 4226 (1995). 

[24] M. Takigawa et al. Phys. Rev. B 55, 14129 (1997). 

[25] M. Matsumura et al. Pliys. Rev. B 56, 8938 (1997). 



References 

[I] N.W. Ashcrott and N.D. Mermin, Solid State Physics, 
Saunders College/Holt, Rinchart and Winston, 1976. 

[2] S. Eggert, I. Affleck and M. Takahashi, Phys. Rev. Lett 73, 
332 (1994). 

[3] A divergent slope as T — > is surprising, because it 
appears to contradict the third law of thermodynamics 
lim^^oAS -* 0. Therefore, |^ = _^|£ = |i| 
should vanish as T — * 0. However, for the spin-1/2 chain 
the derivatives in respect to field and temperature do not 
commute in this limit. 

[4] H.G. Evertz, G. Lana and M. Marcu, Phys. Rev. Lett. 70, 
875 (1993) 

[5] H.G. Evertz, preprint |cond-mat/9707221| (1997). 

[6] U. Wolff, Phys. Rev. Lett. 62, 361 (1989) 

[7] B.B. Beard and U.J. Wiese, Phys. Rev. Lett. 77, 5130 
(1996) 

[8] J.D. Reger, A. P. Young, Phys. Rev. B 37, R5978 (1988). 

[9] A.W. Sandvik, E. Dagotto, D.,I. Scalapino, Phys. Rev. B 
56, 11701 (1997) 

[10] S. Sachdcv, C. Buragohain, M. Vojta, Science 286, 2479 
(1999) 

[II] N. Nagaosa, Y. Hatsugai, M. Imada, J. Phys. Soc. Jpn. 
58, 978 (1989) 

[12] K.H. Hoglund, A.W. Sandvik, Phys. Rev. Lett. 91, 77204 
(2003) 

[13] K.H. Hoglund, A.W. Sandvik, Phys. Rev. B 70, 24406 
(2004) 



